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ABSTRACT 


An  iterative  beam  processor  is  developed  which  in 
the  lim-t  is  identical  to  the  mixed-signal  processor 
(Dean  et.  al.,  19683.  Assuming  that  two  events  arrive 
simultaneously  at  an  array  consisting  of  N sensors, 
the  array  is  first  beamed  on  one  of  the  two  epicenters 
to  produce  a signal  estimate  for  this  event  (O'th 
iteration).  This  signal  estimate  is  then  time-shifted 
and  subtracted  from  each  of  the  original  N seismograms 
in  an  attempt  to  remove  the  signal  from  the  original 
seismograms.  The  new  set  of  records,  each  containing 
N stripped  seismograms,  is  then  beamed  to  produce  a 
signal  estimate  for  the  second  event.  The  signal 
estimate  for  the  second  event  is  now  time-shifted, 
subtracted  from  the  original  N seismograms,  and  the 
stripped  seismogram  are  rebeamed  on  the  first  event. 

The  process  is  repeated  until  differences  in  successive 
signal  estimates  for  the  desired  event  fall  below  a 
predetermined  threshold.  The  iterative-beam  processor 
has  great  practical  (and  intuitive)  appeal.  For  seven 
or  more  elements,  the  iterative  process  converges  in 
a te.  iterations  requiring  only  a few  shift  and  sum 
operation  per  data  point,  while  the  equivalent  mixed- 
signal  (asymptotic  maximum-likelihood)  processor 
requires  a convolution  for  each  data  point. 
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j INTRODUCTION 

{ The  classical  method  for  separating  two  signals 

which  arrive  simultaneously  at  a seismic  array  is  to 
beam  the  array  on  each  of  the  event  epicenters.  Un- 
fortunately, the  simple  beam  does  not  always  yield 
satisfactory  signal  estimates  due  to  contamination  of 
one  signal»s  estimate  by  leakage  from  the  other  signal. 
Shumway  (1972),  however,  demonstrated  that  a mixed- 
signal  (asymptotic  maximum- likel ihood)  processor  yields 
better  signal  estimates  than  does  simple  beamforming. 
Using  various  TFO  subarrays,  and  superimposed  signals 
from  an  earthquake  in  the  Fox  Islands  and  one  in  the 
Tonga  Islands,  Shumway  found  that  the  superiority  of 
the  mixed-signal  processor  was  especially  pronounced 
for  small  arrays.  Similarly,  Cohen  (1972),  in  a study 
of  the  coda  suppression  capabilities  of  the  beam  and 
mixed-signal  processors  performed  using  signals  recorded 
at  TFO,  found  that  the  coda  attenuation  obtained  using 
a 7-element  subarray  and  the  mixed-signal  processor 
was  comparable  to  that  obtained  using  a 19-element 
subarray  and  the  beam. 

While  the  superiority  of  the  mixed-signal  proces- 
sor has  been  established  by  these  and  other  studies, 
it  is  not  always  practical  to  use  the  maximum-likeli- 
hood approach.  A convolution  is  required  for  each 
data  point;  as  such,  when  data  from  a large  array  with 
long  moveouts  must  be  processed,  limitations  on  core 
and  computational- time  may  prevent  its  use.  The 
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iterative-beam  approximation  to  the  mixed-signal  pro- 
cessor can  overcome  these  drawbacks.  Also,  the  itera- 
tive beamforming  approach  is  similar  in  concept  to 
simple  beamforming  and  (by  minimal  investments  in 
additional  system  programming)  may  be  able  to  take 
advantage  of  existing  hardware  and  software  in  opera- 
tional systems. 

In  this  report  we  first  introduce  the  iterative- 
beam  processor,  and  demonstrate  both  theoretically  and 
empirically  that  in  the  limit  it  yields  results  identi- 
cal to  those  produced  by  the  mixed-signal  processor. 

We  next  determine  the  characteristics  of  the  iterative- 
beam  approximation,  specifically  inquiring  into  the 
rate  of  convergence.  Finally,  we  examine  the  action 
of  the  iterative-beam  processor  on  noise.  Here,  the 
data  suggest  that  the  iterative-beam  processor  (and 
hence  the  mixed-signal  processor)  behaves  in  a manner 
similar  to  that  of  the  simple  beam;  that  is,  it  reduces 
the  rms  levels  of  the  final  noise  estimates  by  approxi- 
mately N^,  where  N is  the  number  of  sensors  in  the 
array . 


SIGNAL  ANALYSIS  THEORY 


Consider  the  collection  of  signals 


yjt) 


I ? 


u*-“  k=l 


s,  (u)  + n.rt) 


where 


j “ l***«*n,  is  the  collection  of  N observed 
time  series;  Sj^(t),  k = l,...,p,  is  the  collection  of 
signals  to  be  estimated;  x^j^(t)  = <5(t-T.j^)  [where  6(t)  = 1 
for  t = 0,  zero  otherwise];  and  n. (t)  is  the  noise 
on  the  j'th  channel. 

iVe  assume  that  the  signals  can  be  estimated  by 
linear  estimates  of  the  type: 

■ jli  L (2) 


where  (t)  is  a p x n matrix  of  filter  functions  to 
be  determined. 

Using  this  formulation.  Dean  et.  al.,  (1968)  showed 
that  in  the  frequency  domain,  and  for  uncorrelated 
noise,  the  matrix  of  filter  coefficients  which  produces 
the  best  linear  unbiased  estimates  of  the  Sj^(t)  is 
given  by 

'-1(0))  * [X*(o))  X(o))]’^X*(o)). 
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Transforming  the  M(o))  calculated  at  each  frequency 

back  into  time,  the  signal  estimates  s^(t)  can  then  be 

obtained  by  convolving  the  filters  with  the  y.(t).  For 

a single  signal  (k=’)  in  uncorrelated  noise,  L(t)  is 
the  beam.  ^ 


A block  diagram  of  the  iterative-beam  processor 
is  given  in  Figure  1.  That  this  processor,  in  the 
limit  as  the  number  of  iterations  becomes  large,  yields 
results  identical  to  those  produced  by  the  mixed-signal 
processor  is  shown  in  Appendix  I.  Briefly,  with  respect 
to  Figure  1,  we  assume  that  signals  from  two  events  are 
recorded  at  an  array  consisting  of  N sensors.  The  array 
IS  first  beamed  on  one  of  the  two  epicenters  to  produce 
a signal  estimate  for  the  chosen  event.  This  estimate 
IS  then  time-shifted  and  subtracted  from  each  seismo- 
gram of  the  original  N seismogram  in  an  attempt  to 
remove  the  given  event's  signal  from  the  original  seis- 
mograms. We  now  have  a new  set  of  array  seismograms 
which  contains  primarily  those  signals  corresponding 
to  the  second  signal.  This  set,  containing  N stripped 
seismograms,  is  now  beamed  to  produce  a signal  esti- 
mate  for  the  second  event.  The  subtraction  and  beaming 
process  IS  repeated  until,  for  example,  the  differences 
in  successive  signal  estimates  for  the  desired  event 
falls  below  a pre-determined  threshold. 


The  technique  of  stripping  seismograms  to 
various  arrivals  is  identical  to  the  method  of 
tive  subtraction  of  coherent  noise  described  by 
Passechnik  (1972).  Further,  it  is  of  interest 


enhance 

consecu 

to  note 
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Figure  1.  Iterative-be 


am  processor. 


thac  the  ''stripping”  process  discussed  by  Smart  (1972) 
can  be  seen  to  be  the  first  step  in  this  iterative 
process,  and  thus  to  be  a first  approximation  to  a 
maximum- likelihood  process. 

There  are  actually  two  iterative  methods  for  ob- 
taining a signal  estimate  for  a selected  event.  We 
can  proceed  as  discussed  previously,  beaming  first  on 
the  event  for  which  we  desire  a signal  estimate,  and 
iterating  according  to  the  scheme  shown  in  Figure  1. 
Alternatively,  we  can  start  by  beaming  on  the  second 
signal  present  (for  which  we  may  not  desire  a signal 
estimate),  iterate  according  to  the  method  of  Figure 
1,  but  take  the  signal  estimates  for  the  event  of 
interest  at  the  middle  of  each  iteration.  It  should 
be  noted,  however,  that  the  estimate  of  the  signal 
extracted  at  the  middle  of  an  iteration  is  biased 
(see  Appendix  II  for  details). 

In  our  analysis  of  the  iterative  beam  processor 
corresponding  signal  amplitudes  on  each  channel  of  the 
original  seismograms  are  approximately  equal.  As  such, 
we  weight  the  beam  estimates  by  the  reciprocal  of  the 
number  of  channels  summed.  The  iterative  solution 
computed  using  equal  weighting  factors  corresponds 
directly  to  the  solution  which  would  be  obtained  using 
the  mixed-signal  processor,  since  the  signal  model  in 
these  cases  is  given  by: 

y-jCt)  = s^(t  - T.^)  + S2(t  - T.p  + n.(t) 
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where 


yj(t)»  j ~ is  the  collection  of  N observed 

time  series;  Sj(t)  and  S2(t)  are  propagating  plane 
waves;  and  T.^  are  appropriate  delay  times  for  the 
j ' th  sensor  and  nj(t)  is  the  noise  signal  on  the  j ' th 
sensor. 
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Data 

Seismograms  recorded  at  the  inner  19  elements  of 
TIO  (figure  2,  Table  I)  for  two  events  (Table  II),  one 
from  the  Fox  Islands  (Figure  3),  and  the  other  from 
the  Tonga  Islands  (Figure  4),  were  used  to  investigate 
the  characteristics  and  capabilities  of  the  iterative- 
beam  processor.  Both  events  were  recorded  at  each 
sensor  with  signal-to-noise  ratios  of  25  to  30  db . 

To  simulate  the  near-simultaneous  arrival  of  two 
events  at  TFO,  seismograms  for  those  events  wore 
superimposed;  delays  for  the  arrivals  from  the  Fox 
Islands  event  relative  to  the  arrivals  from  the  Tonga 
Islands  event  were  taken  to  be  on  the  order  of  4 to  5 
seconds.  Thus  the  composite  data  sot  (Figure  5) 
corresponds  closely  to  the  data  used  by  Shumway  (1972) 
to  evaluate  the  mixed-signal  processor. 

I^monstrations  of  the  Mixed-Signal  and  Iterative-Ream 
I rocosso rs  — — 


To  demonstrate  initially  that  the  mixed-signal 
processor  and  iterative-beam  processor  with  equal 
weights  yield  identical  results,  the  processors  wore 
applied  to  the  mixed-signal  data  shown  in  channels 
1-7  of  Figure  5.  The  results  of  the  analyses  are 
shown  in  Figure  6.  As  soon  in  Figures  6b  and  d,  the 
signal  estimates  for  both  events  obtained  using  the 
iterative-beam  processor  are  virtually  identical  to 
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TABLE  I 

Coordinates  of  TFO  Instruments 
Z1  - Z19 
SHORT  PERIOD 


LATITUDE  LONGITUDE  ELEVATION  (Meters) 


Z-1 

34 

16 

42.3300N 

111 

18 

12.2560W 

1615.42 

1-1 

34 

19 

14„1190N 

111 

19 

26.6720W 

1489 .75 

1-1 

34 

18 

42.5760N 

111 

16 

56.0730W 

1515.3 

Z-4 

34 

16 

13.8910N 

111 

14 

56.7380W 

1474.6 

Z-5 

34 

14 

10.2220N 

111 

17 

2.2670W 

1491.9 

Z-6 

34 

14 

55.6970N 

111 

20 

7.0030W 

1509.4 

1-1 

34 

17 

9.1830N 

111 

21 

24.5490W 

1403.2 

Z-8 

34 

21 

42.2900N 

111 

20 

23.5360W 

1805.4 

1-9 

34 

21 

10.0020N 

111 

17 

12.9640W 

1569.4 

1-10 

34 

20 

52.2090N 

111 

14 

4.9130W 

1658.5 

Z-11 

34 

18 

4.8650N 

111 

12 

23.8630W 

1904.9 

Z-12 

34 

15 

49.5070N 

111 

11 

44.7380W 

1528.29 

Z-13 

34 

13 

48.4020N 

111 

13 

48.3160W 

1513.74 

Z-14 

34 

11 

41.4430N 

111 

16 

30.2940W 

1534.0 

Z-15 

34 

12 

8.0760N 

111 

19 

8.4440W 

1487.4 

Z-16 

34 

12 

32.0810N 

111 

22 

13.9550W 

1462.9 

Z-17 

34 

15 

3.9540N 

111 

23 

28.2740W 

1426.4 

Z-18 

34 

17 

40.8760N 

111 

24 

39.9550W 

1664.9 

Z-19 

34 

19 

39.8540N 

111 

22 

32.4010W 

1588.55 
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Figure  3.  Seismograms  for  the  Fox  Islands  Event. 
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Figure  5.  Superposition  of  the  Tonga  and  Fox  Islands  Hvents. 


A 


the  results  produced  by  the  mixed-signal  processor. 
Further,  note  the  significant  reduction  in  leakage 
signals*  on  these  estimates  over  the  results  obtained 
using  simple  beam  sums  (Figures  6 a and  c) . For 
example,  the  leakage  signal  from  the  Tonga  Islands 
event,  which  can  be  observed  prior  to  the  arrival  of 
the  Fox  Islands  signal  in  Figure  6c,  is  reduced  by 
over  8 db  in  the  estimates  shown  in  Figure  6d. 

The  results  from  similar  analyses  of  all  19  TFO 
channels  shov-n  in  Figure  5 aie  shown  in  Figure  7. 

Again,  me  iterative-beam  and  mixed-signal  processors 
yield  practically  identical  results.  Further,  while 
the  leakage  signals  in  the  19-element  beams  (Figures 
7a  and  c)  are  significantly  reduced  over  those  observed 
in  the  7-element  beams  (Figure  6a  and  c) , we  note  that 
the  suppression  of  leakage  signals  using  19-elements 
and  the  simple  beam  is  comparable  to  the  attenuation 
obtained  using  the  7-element  subarray  and  the  iterative- 
beam  or  mixed-signal  processor, 

A further  demonstration  that  the  iterative  approxi- 
mation yields  results  identical  to  .hose  produced  by 
the  mixed-signal  processor  is  shown  in  Figure  8.  Here, 
while  only  the  signals  from  the  Fox  Islands  event 
(Figure  3;  19  elements  analyzed)  are  actually  present, 
we  sought  estimates  of  signals  from  both  the  Fox  and 
Tonga  Islands.  The  "estimates”  for  the  Tonga  Islands 

signal  is  that  proportion  of  Signal  1 
which  leaks  into  our  estimate  for  Signal  2 (and  vise 
versa) . 
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TOP:  ITERATIVE  ESTIMATE  FOR  FOX  ISLANOS  SIGNAL(4  ITERATIONS) 
BOTTOM:  MIXEO-SIGNAL  ESTIMATE  FOR  FOX  ISLANOS  SIGNAL 


signal  are,  of  course,  leakage  signals.  As  seen  in 
Figures  8b  and  d,  the  estimates  obtained  using  the 
iterative-beam  and  mixed-signal  processors  are  vir- 
tually identical.  Further,  the  simple  beam  estimate 
for  the  Fox  Islands  signal  (Figure  8c3  appears 
identical  to  the  estimates  obtained  using  the  more 
sophisticated  processors. 

A more  quantitative  analysis  of  the  leakage  sig- 
nals (Figures  8a  and  b)  can  be  made  by  comparing  the 
power  in  these  signals.  We  justify  this  approach  on 
the  grounds  that  the  leakage  represents  an  additional 
noise  signal  which  would  appear  in  the  estimate  of 
the  Tonga  Islands  signal,  had  one  been  present.  Using 
a 35-second  window  beginning  at  a point  corresponding 
to  the  first  arrival  of  the  Tonga  Islands  signal 
(Figures  8c  and  dj  , a measure  of  the  power  IT  in  the 
various  estimates  of  the  leakage  signal  is  given  as 
follows : 

K = [T-l  J s*2(t)itj 

t-1 

A 

where  s(tj  is  the  estimate  for  the  leakage  signal,  and 
T is  the  length  of  the  sample  analyzed. 

For  the  simple  beam  (Figure  8a),  ^ is  0.914mu^, 
while  for  the^iterative-beam  and  mixed-signal  estimates 
(Figure  8b),  R is  found  to  be  0.174  and  0.176mu^, 
respectively.  Thus,  in  this  case,  the  more  sophisti- 
cated processors  yield  an  additional  7.2  db  of  coda 
attenuation  over  that  provided  by  the  simple  beam.  In 
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general,  the  coda  attenuation  capability  of  the  mixed- 
signal  processor  (and  its  approximation,  the  iterative- 
beam  processor)  over  that  obtained  using  the  simple 
beam  is  on  the  order  of  3 to  5 db,  although  signifi- 
cantly more  improvement  can  be  expected  if  the  signals 
are  close  in  velocity  space  (Cohen,  1972). 

Characteristics  of  the  Iterative-Beam  Processor 

Having  demonstrated  that  the  iterative-beam  pro- 
cessor yields  solutions  identical  to  those  produced 
by  the  mixed-signal  processor,  we  now  examine  some 
characteristics  of  the  iterative  approximation.  Spe- 
cifically, we  discuss  rate  of  convergence  as  a function 
of  how  the  processor  is  applied,  and  the  number  of 
channels  analyzed  (i.e,,  the  array  size). 

Two  Signals  - 7 Channels 

Application  of  the  iterative-beam  processor  to 
the  superimposed  Tonga  and  Fox  Islands  signals  recorded 
on  the  seven  inner  elements  of  TFO  (Figure  5,  Channels 
1-7)  yields  the  results  given  in  Table  III  and  shown 
in  Figures  9 through  12  (all  signal  estimates  shown 
are  normalized  to  the  same  peak  amplitude).  Figures 
9a  and  11a  suggest  that  the  two  iterative  processes 
are  converging  to  the  same  answer,  but  that  one,  for 
which  the  initial  beam  is  on  the  event  of  interest,  is 
converging  more  quickly  than  the  other.  Even  removing 
the  predictable  bias  on  the  mid- iteration  signal  esti- 
mates (Figures  9b  and  11b)  does  not  completely  correct 


-20- 


f 


TABLE  III 

Iterative-Beam  Analysis 
Two-Signals  7-Channels 


Signal  1 (Tonga  Islands) 

Signal  2 (Fox 

Islands) 

Initial  Beam  on  Tonga  Is. 

Initial  Beam 

on  Fox  Is. 

Maximum  Amp. 

Iteration  Sum  Trace  (my) 

Iteration 

Maximum  Amp. 
Sum  Trace  (my) 

0 13.12 

0 

14.54 

1 12.93 

1 

14.69 

2 12.90 

2 

14.79 

3 12.90 

3 

14.89 

4 12.91 

4 

14.90 

Initial  Beam 

on  Fox  Is. 

Initial  Beam 

on  Tonga  Is. 

Iteration 

Maximum  Amp 
Sum  Trace  (my) 

Iteration 

Maximum  Amp. 
Sum  Trace  (mu) 

1 

9.48 

1 

11.35 

2 

11.31 

2 

13.61 

3 

12.00 

3 

14.32 

4 

12.33 

4 

14.62 

5 

12.52 

5 

14.73 
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ITERATION 

^iguro  9a.  Convergence  characteristics  for  the  Tonga  Island 
Iterative-beam  signal  estimates,  superimposed  Tonga  and  Fox 
Islands  signals,  7 channels. 
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ITERATION 


Figure  9b.  Convergence  characteristics  for  the  Tonga  Islands 
Iterative-beam  signal  estimates,  superimposed  Tonga  and  Fox 
Islands  signals,  7 channels,  predictable  bias  removed. 
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figure  11a.  Convergence  characteristics  for  the  Fox  Islands 
iterative-beam  signal  estimates,  superimposed  Tonga  and  Fox 
Islands  signals,  7 channels. 
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ITERATION 


Figure  11b.  Convergence  characteristics  for  the  Fox  Islands 
iterative-beam  signal  estimates,  Tonga  and  Fox  Islands 
signals,  predictable  bias  removed. 


(SIMPLE  BEAM) 


for  the  contamination  of  a given  signal  by  its  own 
echoes.  Further,  as  seen  in  Figure  10b,  the  mid- 
iteration signal  estimates  for  the  Tonga  Islands  event 
(initial  beam  on  the  Fox  Islands)  erl.ibit  precursors 
to  the  first  arrival  which  derive  from  echoes  of  the 
‘ Tonga  Islands  signal.  It  should  be  noted  that  the 

large  precursors  on  the  signal  estimation  for  the  Fox 
Islands  event  (Figure  12a)  derive  from  leakage  of  the 
Tonga  Islands  signal,  and  not  from  echoes  of  the  Fox 
Islands  event  (though  these  must  be  present). 

I The  results  obtained  here  amplify  our  earlier 

I comments  on  the  proper  mode  for  implementation  of  the 

[ iterative-beam  processor.  That  is,  one  should  beam 

initially  on  the  signal  of  interest,  and  use  only 
those  signal  estimates  produced  at  the  end  of  a com- 
plete iteration. 

Choice  of  a test  for  convergence  of  the  iterative 
process  is  somewhat  subjective.  For  many  qualitative 
purposes  (e.g.  detections),  the  first  iteration  is 
significantly  improved  over  the  beam  to  such  an  extent 
that  the  process  can  be  terminated  without  a test.  In 
general,  however,  a quantitative  decision  has  to  be  made 
as  to  what  constitutes  a "significant  improvement". 
Though  the  analyses  discussed  in  this  report  were  arbi- 
trarily terminated  on  the  4^^  or  5^^  iteration,  we 
could  have,  for  example,  computed  the  rms  amplitude 
of  the  difference  between  successive  signal  estimates, 
and  terminated  the  process  when  the  rms  value  fell 
below  a prescribed  level. 
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For  the  7-channel,  iterative-beam  computations 
discussed  here,  the  total  CPU  time  for  the  Tonga  and 
Fox  Island  signal  estimates  (four  iterations  each)  on 
an  IBM  360/44  was  on  the  order  of  13  minutes.  For 
both  signals  and  only  one  iteration  on  each  solution, 
the  time  is  approximately  3.0  minutes.  This  is  one- 
half  the  CPU  time  required  for  the  corresponding  solu- 
tions using  the  mixed-signal  processor  (-6  minutes). 

One  must  also  remember  that  for  more  than  7 elements, 
the  iterative  process  will  converge  even  more  accu- 
rately in  two  iterations  for  any  number  of  channels, 
while  the  number  of  convolutions  needed  per  point  for 
the  maximum  likelihood  approach  increases  in  propor- 
tion to  N,  the  number  of  channels.  As  a result,  the 
maximum- likel ihood  processor  is  impractical  for  LASA 
short-period  data,  while  iterative  beamforming  is 
practical. 

Two  Signals  - 19  Channels 

For  the  superimposed  Tonga  and  Fox  Islands  signals 
recorded  on  a 19-element  TFO  array  (Figure  5)  , the 
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iterative-beam  processor  yields  the  results  given 
Table  IV  and  shown  in  1‘igures  13  through  16. 

The  total  CPU  time  for  the  signal  estimates  ( 
iterations  on  each  signal)  was  on  the  order  of  30 
However,  because  the  iterative  process  converged  s 
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TABLE  IV 

I terat ive- Beam  Analysis 
Two-Signal  19-Channels 


Signal  1 (Tonga  Islands) 
Initial  Beam  on  Tonga  Is. 


I teration 

0 

1 

2 

3 

4 


Maximum  Amp 
Sum  Trace  (mp) 

12.18 

12.48 

12.53 

12.54 
12.53 


Signal  2 (Fox  Islands) 
Initial  Beam  on  Fox  Is. 


I teration 

0 

1 

2 

3 

4 


Maximum  Amp. 
Sum  Trace  (mp) 

14.21 

14.50 

14.50 

14.52 

14.; 


Initial  Beam  on  Fox  Is. 


I teration 


1 

2 

3 

4 

5 


Maximum  Amp. 
Sum  Trace  (mp) 

11.30 

12.18 

12.38 

12.44 

12.47 


Initial  Beam  on  Tonga  Is. 
I teration 


1 

2 

3 

4 

5 


Maximum  Amp. 
Sum  Trace  (mp 


13,59 

14,35 

14.52 

14.56 

14.56 
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BEAMED  FIRST  ON  TONGA  ISLANDS 


BEAMED  FIRST  ON  FOX  ISLANDS 
(MID  ITERATION  RESULTS) 


ITERATION 


Figure  13.  Convergence  characteristics  for  the  Tonga  Islands 
iterative-beam  signal  estimates,  superimposed  Tonga  and  Fox 
Islands  signals,  19  channels. 
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BEAMED  FIRST  ON  FOX  ISLANDS 


BEAMED  FIRST  ON  TONGA  ISLANDS 
(MID  ITERATIONIRESULTS) 


ITERATION 


Figure  15.  Convergence  characteristics  for  the  rox  Islands 
iterative-beam  signal  estimates,  superimposed  Tonga  and  Fox 
Islands  signals,  19  channels. 
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at  7 minutes.  This  compares  favorably  with  the  CPU 
time  required  for  corresponding  solutions  using  the 
mixed-signal  processor  (~12  minutes  for  19  channels^ . 

One  Signal  - 7 Channels 

Table  V and  Figures  17  through  20  show  the  results 
of  processing  the  Fox  Islands  signal  (Figure  33  under 
the  assumption  that  two  signals  - one  from  tlic  Fox 
Islands  and  one  from  the  Tonga  Islands  - are  present. 

A plot  of  the  sum  trace  rms  amplitude  for  the  non- 
existent Tonga  Islands  signal  is  shown  in  Figure  17 
as  a measure  of  the  portion  of  the  Fox  Island  events 
which  would  leak  into  our  estimates  of  the  Tonga 
Islands  signal  (had  a Tonga  Island  signal  been  present^. 

In  this  case,  beaming  first  on  the  Fox  Islands 
and  extracting  the  Tonga  Islands  signal  (noise) 
estimates  at  mid- iteration  yields  solutions  for  the 
Tonga  Islands  event  which  stabilize  more  quickly  than 
those  produced  by  bearing  first  on  this  non-existent 
event.  This  is  not  surprising  when  one  considers  that 
the  estimates  for  the  non-existent  Tonga  Islands  event 
which  are  stripped  out  when  one  beams  first  on  this 
event  consist  of  a combination  of  noise  and  leakage 
from  the  Fox  Islands  event.  By  beaming  first  on  the 
Fox  Islands  event  and  stripping  out  its  signal, 
leakage  into  the  Tonga  Islands  estimate  is  minimized 
in  the  early  iterations. 

These  data  suggest  that  if  it  is  not  known  whether 
a second  signal  is  present,  or  under  low  signal- to-n'>ise 
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TABLE  V 

Iterative-Beam  Analysis 
One-Signal  7-Channels 


Signal  1 (Tonga  Islands) 
Initial  Beam  on  Tonga  Is. 

Sum  Trace 

Iteration  rms  (my) 


Signal  2 (Fox  Islands) 
Initial  Beam  on  Fox  Is. 

Maximum  Amp. 

Iteration  Sum  Trace  (my) 


0 

1 

2 

3 

4 


2.311 

1.108 

0.782 

0.697 

0.676 


0 

1 

2 

3 

4 


IS. 34 
15.20 
15.17 
15.15 
15.19 


Initial  Beam  on  Fox  Is. 

Sum  Trace 

Iteration  rum  (my) 


Initial  Beam  on  Tonga  Is. 

Maximum  Amp. 

Iteration  Sum  Trace  (my) 


1 

2 

3 

4 

5 


0.494 

0.591 

0.628 

0.646 

0.657 


1 

2 

3 

4 

5 


11.94 

14.16 

14.77 

14.96 

15.06 
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ITERATION 


ligurc  17.  Convergence  characteristics  for  the  Tonga  Islands 
iterat  i ve-lieam  signal  estimates,  I'ox  Islands  signal,  7 channels. 
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Figure  18,  Iterative-beam  Tonga  Islands  signal  estimates 
Fox  Islands  signal,  7 channels. 
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Figure  19.  Convergence  characteristics  for  the  Fox  Islands 
iterative-beam  signal  estimates.  Fox  Islands  signal,  7 channels. 
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conditions  for  the  second  event  (second  events  signal 
buried  in  the  coda  of  a large  event) , one  should  beam 
first  on  the  known  event's  epicenter,  and  examine  the 
i icrat ion  signal  estimates  for  the  second  event. 

If  a weak,  second  event  is  indeed  present,  its  signal 
estimate  will  be  biased,  but  its  waveform  will  stabi- 
lize more  rapidly  than  if  we  beam  first  on  this  signal. 

Thus,  when  the  mere  existence  of  a weak  second  signal, 
rather  than  its  amplitude  and  detailed  waveform,  is  of 
interest,  it  would  seem  best  to  initially  beam  on  the 
first  signal. 

One  Signal  - 19  Channels 

Application  of  the  iterative-beam  processor  of 
the  Fox  Islands  signal  recorded  on  a 19-element  TFO 
array  (Figure  3)  yields  the  results  shown  in  Table  VI 
and  Figures  21  through  24,  Except  for  a scale  factor, 
the  various  signal  estimates  for  the  Fox  Islands  event 
are  almost  indistinguishable  from  one  another.  Due  to 
the  large  number  of  array  elements  (and  hence,  in  this 
case,  the  large  aperture  of  the  array)  signal  estimates 
for  the  non-existent  Tonga  Islands  event  stabilized 
after  one  iteration.  If  a weak  Tonga  Islands  signal 
had  been  present,  however,  it  would  be  best  for 
detection  purposes  to  extract  the  mid- iteration  signal 
estimates  for  this  event  produced  by  beaming  first  on 
t’’e  Fox  Islands  event. 
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TABLE  VI 

Iterative-Beam  Analysis 
One-Signal  19-Channels 


Signal  1 (Tonga  Islands) 

Signal  2 (Fox 

: Islands) 

Initial  Beam  on  Tonga  Is. 

Initial  Beam 

on  Fox  Is. 

Sum  Trace 

Iteration  rms  (mg) 

I teration 

Maximum  Amp. 
Sum  Trace  (my) 

0 0.956 

0 

14.27 

1 0.447 

1 

14.29 

2 0.423 

2 

14.30 

3 0.421 

3 

14.29 

4 0.421 

4 

14.26 

Initial  Beam 

on  Fox  is. 

Initial  Beam 

on  Tonga  Is. 

Iteration 

Sum  Trace 
rms  (mg) 

Iteration 

Maximum  Amp. 
Sum  Trace  (mg) 

1 

0.385 

1 

13.81 

2 

0.405 

2 

14.11 

3 

0.412 

3 

14.18 

4 

0.415 

4 

14.21 

5 

0.418 

5 

14.24 
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B(AM(D  riRSI  ON  lUNCA  ISIANOS 


BEAMIO  FIRST  ON  FOX  ISLANDS 
IMIO  ITFRATION  RFSUlTSl 


2 3 

ITERATION 


ij;urc  21.  Convergence  cliaracterist  ics  for  the  longa  Islands 
i tc  rat  i ve-beum  signal  estimates,  fox  Islands  signal,  ID  cliannc 
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Figure  22.  Iterative-beam  Tonga  Islands  signal  estimtes,  Fox 
Islands  signal,  19  channels. 


ITERATION 


Hgurc  23.  Convergence  characteristics  far  the  !'ox  Islands 
iterative-beam  signal  estimates,  I ox  Islands  signal,  19  channels 


I 
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ITERATION 


Figure  24,  Iterative-beam  Fox  Islands  signal  estimates,  Fox 
Islands  signal,  19  channels. 


Noise 


Let  us  now  examine  the  signal  estimates  obtained 
using  the  beam,  iterative-beam,  and  mixed-signal  pro- 
cessors in  cases  where  only  noise  is  present.  The 
first  noise  sample  analyzed  was  recorded  prior  to  the 
arrival  of  the  Fox  Islands  signal.  Application  of  the 
processors  to  the  noise  recorded  on  7 and  19  channels 
yields  the  results  listed  in  Tables  VII  a and  b,  and 
shown  in  Figures  25  and  26,  respectively.  In  every 
case,  the  signal  estimates  pro«-aced  on  the  fourth 
ite.*ation,  regardless  of  the  event  beamed  initially, 
were  identical.  The  noise  waveforms  converge  as  quickly 
as  the  signal  waveforms  that  would  be  expected  from  the 
linearity  of  the  process. 

Analyses  of  the  noise  recorded  prior  to  the  arri- 
val of  the  Tonga  Islands  signal  yield  similar  results. 
For  that  reason,  only  the  tabulated  values  for  the  beam 
and  iterative-beam  computations  are  given  here  (Tables 
VIII  a and  b) . 

Of  primary  concern  to  us  is  a comparison  of  the 
rms  levels  for  the  beam  and  iterative-beam  noise 
estimates.  We  might  expect  the  iterative-beam  pro- 
cessor to  be  similar  to  the  simple  beam  when  processing 
noise;  that  is,  we  expect  both  processors  to  reduce 
the  rms  levels  of  the  final  noise  estimates  by  some- 
thing on  the  order  of  the  N^,  where  N is  the  number 
of  channels  processed. 

Table  IX  shows  a comparison  of  the  rms  levels  for 
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TABLE  Vila 

Iterative-Beam  Analysis 
Noise  (Fox  Islands)  7-Channels 


Signal  1 

(Tonga  Islands) 

Signal  2 (Fox 

Islands) 

Initial  Beam  on  Tonga  Is. 

Initial  Beam 

on  Fox  Is. 

Iteration 

Sum  Trace 
rms  (mu) 

Iteration 

Sum  Trace 
rms  (my) 

0 

0.282 

0 

0.290 

1 

0.250 

1 

0.274 

2 

0.236 

2 

0.273 

3 

0.230 

3 

0.274 

4 

0.229 

4 

0.277 

Initial  Beam 

on  Fox  Is. 

Initial  R ;am 

on  Tonga  Is 

Iteration 

Sum  Trace 
rms  (my) 

Iteration 

Sum  Trace 
rms  (my) 

1 

0.101 

1 

0,103 

2 

0.144 

2 

0.149 

3 

0.170 

3 

0.181 

4 

0.186 

4 

0.206 

5 

0.198 

5 

0.225 
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TABLE  VI Ib 

Iterative-Beam  Analysis 
Noise  (Fox  Islands)  19-Channels 


I 

I 


Signal  1 (Tonga  Islands) 

Signal  2 (Fox 

Islands) 

Initial  Beam 

on  Tonga  Is. 

Initial  Beam 

on  Fox  Is. 

Iteration 

Sum  Trace 
rms  (mu) 

I teration 

Sum  Trace 
rms  (mu) 

0 

0.187 

0 

0.204 

1 

0.158 

1 

0.192 

2 

0.150 

2 

0.189 

3 

0.147 

3 

0.188 

4 

0.147 

4 

0.187 

Initial  Beam 

on  Fox  Is. 

Initial  Beam 

on  Tonga  Is 

Iteration 

Sum  Trace 
rms  (mu) 

Iteration 

Sum  Trace 
rms  (mu) 

1 

0.101 

1 

0.116 

2 

0.121 

2 

0.150 

3 

0.129 

3 

0.166 

4 

0.134 

4 

0.174 

5 

0.138 

5 

0.179 
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(f)  MIXED-SIGNAL  PROCESSOR;  FOX  ISLANDS 


Figure  26,  Conparison  of  noise  estimates  from  noise  recorded 
prior  to  the  arrival  of  the  Fox  Islands  signal,  19  channels. 


TABLE  Villa 

Iterative-Beam  Analysis 
Noise  (Tonga  Islands)  7-Channels 


Signal  1 (Tonga  Islands)  Signal  2 (Fox  Islands) 


Initial  Beam 

on  Tonga  Is. 

Initial  Beam 

on  Fox  Is. 

1 tcrat ion 

Sum  Trace 
rms  (mu) 

Iteration 

Sum  Trace 
rms  (mp) 

0 

0.216 

0 

0.182 

1 

0.220 

1 

0.157 

2 

0.227 

2 

0.154 

3 

0.232 

3 

0.157 

4 

0.237 

4 

0.162 

Initial  Beam 

on  Fox  Is. 

Initial  Beam 

on  Tonga  Is 

Iteration 

Sum  Trace 
rms  (mp) 

Iteration 

Sum  Trace 
rms  (mp) 

1 

0.115 

1 

0.078 

2 

0.161 

2 

0.106 

3 

0.188 

3 

0.124 

4 

0.205 

4 

0.137 

5 

0.218 

5 

0.146 
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TABLh  VI  11b 

Iterative-Beam  Analysis 


Noise  (Tonga  Islands)  19-Channels 


i 

I 


Signal  1 (Tonga  Islands) 
i.'.Ltial  Beam  on  Tonga  Is. 


Signal  2 (Fox  Islands) 
Initial  Beam  on  Fox  Is. 


Sum  Trace 

Iteration  rms  (mu) 


Sum  Trace 

Iteration  nns  (my) 


0 

0.143 

0 

0 102 

1 

0.151 

1 

0.091 

2 

0.158 

2 

0.097 

3 

0,162 

3 

0.102 

4 

0.165 

4 

0.105 

Initial  Beam 

on  Fox  Is. 

Initial  Beam 

on  Tonga  Is 

Iteration 

Sum  Trace 
rms  (my) 

I teration 

Sum  Trace 
rms  (my) 

1 

O 

O 

1 

O 

O 

2 

0.139 

2 

0.088 

3 

0.153 

3 

0.097 

4 

0.159 

4 

0.102 

5 

0.163 

5 

0.105 
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Comparison  of  rms  Noise  estimates 
Two-Sided  F-Test  (a=<,10) 
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the  noise  estimates  produced  by  the  beam  and  iterative- 
beam  processors.  To  tost  the  significance  of  the 
differences  observed,  wo  use  a two  sided  F-test.  The 

degrees  of  freedom  for  the  noise  ostimatof  is  given 
as  follows : 

degrees  of  freedom  = 2BT 

where  B is  the  bandwidth  of  the  noise,  and  T is  the 
sample  length.  For  B . 1 Hz,  and  a sample  length  of 
35  seconds,  2BT  = 70.  If  we  take  a nominal  value  of 
60  for  the  degrees  of  freedom,  recourse  . the  standard 
tables  for  F(05,v^,V2)  fields  a value  c ..53  for  the 
two-sided  F-test  with  a = 0.10.  Only  one  value  for 
the  rms  ratios,  that  of  the  19  channel  Tonga  Islands 
noise  estimate  computed  using  noise  recorded  prior  to 
the  Fox  Islands  event,  exceeds  this  value  (1.57  versus 
1.53);  the  significance  of  the  difference  in  the  rms 
noise  values  for  this  case  may  be  due  to  coherent 
noise  propogating  across  the  array  from  the  direction 
of  the  Tonga  Islands  region.  However,  fo'  a = 0.10, 
we  would  expect  1 out  of  10  ratio  values  to  exceed  * 
the  F-test  value.  Thus  the  data  of  Table  IX  indicate 
that  the  rms  values  for  the  beam  and  iterative-beam 
noise  estimates  are  not  significantly  different.  Fur- 
ther, as  the  19-element  TFO  array,  and  possibly  the 
7-element  array  as  well,  can  be  expected  to  reduce 
the  rms  noise  level  on  the  beamed  trace  by  a factor  of 
N over  the  rms  noise  level  on  the  individual  channels, 
so  too,  apparently,  does  the  iterative  beam  (and. 
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hence,  the  mixed-signal  processor)  reduce  the  rms  level 
of  the  noise  estimates  by  the  N**.  However,  theoretical 
considerations  suggest  that  this  may  not  be  true  if 
the  two  signals  are  close  together  in  velocity  space. 
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CONCLUSIONS 


In  the  limit,  the  iterative-beam  processor  con- 
verges to  the  mixed-signal  processor.  Further,  the 
iterative-beam  processor  has  great  practical  (and 
intuitive)  appeal.  For  se»^en  or  more  array  elements, 
the  iterative  process  converges  in  a few  iterations, 
requiring  only  a few  shift  and  sum  operations  per 
data  point,  while  the  mixed-signal  (maximum  likelihood) 
processor  requires  a convolution  for  every  data  point. 

For  many  qualitative  purposes,  one  iteration 
yields  satisfactory  convergence. 

When  two  signals  are  present,  each  having  approxi- 
mately the  same  amplitude,  one  should  beam  first  on 
the  signal  of  interest,  and  use  only  those  signal  esti- 
mates produced  at  the  end  of  a complete  iteration.  If 
one  signal  is  considerably  stronger  than  the  other, 
however,  the  initial  beam  should  be  on  the  event 
corresponding  to  the  larger  signal,  and  the  estimate 
for  the  smaller  signal,  although  biased,  extracted  at 
mid- iteration. 

The  data  suggest  that  the  iterative-beam  processor 
(and  hence,  the  mixed-signal  processor)  reduces  the 

rms  noise  level  for  traces  well  separated  in  velocity 

u 

space  by  the  N , where  N is  the  number  of  channels. 
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APPENDIX  I 

Theoretical  Proof  of  the  Equivalence  of  Iterative 
Two“Signal  Beaming  and  Maximum-Likelihood  Processing 
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The  Fourier  Transform  of  the  data  on  the  k'th 
channel  is  given  by 

= Sj(„)  C » S^Cu.)  c '^2  (A-i) 

where  j = /^, 


IS : 


Then  the  first  estimate  of  S given  {Y. 

K K “ 1 ^ 

• 1 ? 


^ I e 

k=l 


YkM 


(A-2) 


N 

I 

k«l 


where  H (w)  = A ^ 


IT 


kl 


The  first  estimate  of  is  defined  in  terms  of 


the  residuals 


sU) 


1 ^ 

k-1 


N 


j-(T 


£«1 


] (A-3) 


N 
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where  = ^ [ 

N j<^('r,.,-T.  _) 

Note  that  A = e as  defined  by’  Shumway 

(1972). 

The  second  estimate  of  is  defined  as  a beam  or 
the  residuals 

(V,w  - s(1)m 

k*i  a=i 


where  = i I e^“’^''*  - A(w)l|(*>  („)  ) . 

If  we  rewrite  II^J'  in  terms  of  we  obtain 

"zk^  ' R I - A*(w)ll{J^J. 

These  recursions  arc  general,  and  thus,  one  may 
write : 

"1’'*^^  = ^ lUi  - 


(A-6) 


-Al-3- 


and 


(2) 


‘•2''^  * W 1“2  ■ 


where 

1 


hCJ^Co)) 


h?J5)  (-) 


ano 


U. 

1 


j»Tii 


Note  that  A * U^. 


Combining  (1)  and  (2),  we  obtain: 


. 1 (Uj  - ^ A(U2  • A*Hp))) 

• ^ [Ui  - ^ AU^l  . 

It  can  be  shown  that  this  sequence  converges 
geometrically  to 

.,(-)  - ''“2) 

1 l-N-^|A|^  • 

Further, 

N’^[U,  - N'^  A*U,] 

— T5 2 — — 

2 1 - N ^\A\^ 


(A-7) 


(A-8) 


(A-9) 
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Equations  (A-8)  and  (A-9)  are  soon  to  bo  tho  two- 
^.Snal  «xi.u..IUoii„ood  voctor  equations  (for  spa- 
tially  white  noise)  ^ 

The  error  estimate  after  k itemfinne  • 

K Iterations  is  given  by: 


whe  re 


c ' ^ Of 


pO<»l)^ 


1-p 


T 


(A-IO) 


Thus,  tho  rato  of  conuorsento  is  a function  of  tho 

eparauon  of  tho  piano  ware  vectors  of  tho  two  signals 

act,  C IS  tho  nornalitod  correlation  between  tho 
two  vectors. 

For  comparison,  the  two-signal  maximum- 1 ikel ihood 
Stsnal  estimates  (for  spatially  white  noise)  are 
given  as  follows  (Shumway,  1972)  ; 


"ik  “ ^ [N  e'^  - A(w) 


(A-11) 


and 


"2k  ‘ (N  o-'“’''<2 


j wT.  , 

- A*((aj)  c 


], 


whore  A(„)  - . |a(„,|2  ^ ^ ^ 

1=1 


" j“fTu-T„) 


Then; 


S^Cw) 

SjM 


N 

k*l 


(uj) 


k*l 


-AI-6- 


APPHNDIX  11 

Discussion  of  Bias  in  an  Alternative  Iterative 
Signal  Processor 
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To  see  that  bias  exists  in  the  signal  estimate 
obtained  in  the  alternative  manner  discussed  in  the 
beginning  of  the  text  of  this  report,  consider  the 
following  mixed-signal  model  for  a two-channel  array: 

y^(t)  * Sj(t-Tj^)  + S2(t-T^2) 
y2(t)  = s^(t-l2^)  + S2(t-T22). 

Let  us  assume  that  we  desire  a signal  estimate  for 
Event  1,  but  that  we  will  beam  first  on  the  epicenter 
for  Event  2: 

s"2(t)  ■ SjU)  * ^[SiCt-Tu*Tj2)  ♦ Sj(t-T2i*T22). 


This  is  the  O'th  iteration  for  signal  2.  Shifting 
and  subtracting  this  estimate  to  eliminate  signal  2 
from  the  original  recordings  (y^Ct)  and  y2(t)),  and 
then  beaming  on  Event  1 yields  the  signal  estimate 

Sj(t)  « Sj(t)  - ^[2Sj(t)  * Sj(t-T2j*T22*Tjj-Tj2) 


Si(t+T2i-T22-Tii+Ti2) 1 


We  see  that  this  estimate  for  signal  1 is  distorted 
only  by  echos  of  itself.  Further,  it  contains  a pre- 
dictable, on-beam  contribution  from  signal  1 which 


-AlI-2- 


reduces  the  scale  of  the  signal  estimate;  thus,  the 
estimate  is  biased.  Here,  this  component  has  a value 
of  -j  Sj(t).  In  general,  we  believe  that  the  pre- 
dictable bias  is  given  by: 

where  N is  the  number  of  channels,  m is  an  odd  number 
(m=2p-l,  where  p is  the  iteration  number),  and  s.(t) 
is  the  signal  estimated  at  the  middle  of  the  iteration. 
This  formula  has  not  been  formally  derived,  but  is 
correct  for  all  cases  we  have  worked  out,  including 
N-2,  m=l,3,5,7;  N=3,4,  m=l,3;  and  m=l,  all  N. 

At  the  end  of  the  first  iteration,  the  signal 
estimate  for  livent  2 again  is  of  the  form: 

ft 

S2(t)  = S2(t)  + f(s^(t)]. 

The  above  results  suggest  that  the  signal  esti- 
mates extracted  at  the  middle  of  an  iteration  will 
converge  slowly.  As  such,  it  would  appear  that  the 
proper  mode  for  implementation  of  the  iterative-beam 
processor  if  one  desires  an  accurate  signal  estimates 
is  to  beam  first  on  the  signal  of  interest,  and  to 

use  only  those  signal  estimates  produced  at  the  end 
of  each  iteration. 


-MI-3- 

! 1 

I i 


